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Abstract. 

The total cross section for Coulomb photoproduction of neutral pion pairs in the 
reaction ^yA — > ir°ir°A is estimated within the effective chiral lagrangian approach. 
The amplitude of 77* — > 7r°7r° with one off-shell photon is calculated at 0(p 6 ) in 
the momentum expansion; in addition, nuclear absorption is taken into account. 
Besides its experimental feasibility, the results of the calculation demonstrate that 
the reaction jA — > 7r°7r°A is a powerful source of information on the process 
77 — > 7r 7r° close to threshold. 
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1. Introduction 

The nucleus Coulomb interaction method is an effective approach to study in high- 
energy experiments the low-energy electromagnetic properties of pions in the process 
77 — > mi. There are two possible ways to realize the Coulomb interaction method 
in experiment: the radiative scattering ttA — > Tc-fA of a pion on a nucleus and 
the photoproduction of pion pairs in the Coulomb field of a nucleus 7A — > ttttA. 
The experiments under discussion are based on the fact that for sufficiently small 
momentum transfers the interaction of high-energy particle with nuclei is extremely 
peripheral and thus dominated by scattering on the virtual photons of the Coulomb 
field of the nucleus. The reaction tt A — > ti^A allows to investigate the 77 — > tttt 
process in the region m% w = (p n i + p^) 2 < and to measure the charged pion 
polarizability [p]] at the point m% w = 0, while the reaction 7A — > txtxA occurs in the 
region m^ n > 4m^. close to threshold. In the case of charged pions the two reactions 
give complementary information about the process 77 — > 7r + 7r~ in both the physical 
and the non-physical region. In the case of neutral pions the photoproduction of 
a neutral pion pair in the Coulomb field of a nucleus provides a new source of the 
information on the process 77 — > 7r°7r°, which was measured in the electron-positron 
experiment by the Crystal Ball Collaboration and is under active discussion in 
the context of the physical programme at DA$NE ||. 

The present theoretical interest in the elementary process 77 — > 7r°7r° is 
caused both by experimental data |2| mentioned above and recent progress in 
Chiral Perturbation Theory (ChPT) M up to and including 0(p e ). The process 
77 — > 7T°7T° is very sensitive to the higher-order contributions in ChPT since the first 
non- vanishing amplitude arises from meson loops at 0(p 4 ) without counterterms. 
However, the one- loop amplitude at 0(p 4 ) calculated in paper || does not describe 
the data even near threshold ||. This is not surprising, as the analysis of the 
77 — > 7r 7T° amplitude based on dispersion relations [|7| demonstrates the importance 
of unitarity corrections corresponding to higher orders, next to 0(p A ). In fact, 
the two- loop calculation at 0{p & ) carried out in paper || gives a considerable 
improvement of the description within ChPT. For completeness we mention that 
there is also another consideration of 77 — > 7r°7r° up to one-loop order corresponding 
to 0(p 5 ) in the treatment of Generalized ChPT ||. 

The contribution from the 0(p 6 ) counterterms was estimated in paper || |f]from 
the low-energy meson phenomenology with resonance exchange saturation. In papers 



10| , |lT| the counterterms were fixed as effective meson lagrangians with higher- 
order derivative terms obtained from the bosonization of Nambu-Jona-Lasinio model 
(NJL). The sensitivity of 77 — > 7r°7r to higher-order corrections including Born 

f The same paper presents a diagrammatical expansion of the amplitude 77 — * tt°tt q up to 0(p 6 ). 



3 



contributions from the effective meson lagrangian at 0(p 6 ) makes this process a 
valuable source of the experimental information essential for the test of bosonized 
chiral lagrangians at 0(p 6 ). 

In this paper we investigate the possibility of studying the process 77 — > 7r°7r° 
near threshold in the photoproduction of neutral pion pairs in the Coulomb field 
of a nucleus. For the first time the reaction jA — > ttttA was considered in this 
context in paper [12] with a one-loop amplitude of 77 — > 7r°7r° at 0(p 4 ). The total 



cross sections 7 A — > txtxA were estimated for energies of the incident photon at 20 
GeV and 40 GeV for momentum transfer cutoffs 5 MeV and 10 MeV, respectively. 
However, due to the large incident energy and small cutoff, the nuclear absorption 
and the offshellness of the virtual photon from the nuclear Coulomb field were not 
taken into account. 

In this note, we extend and improve the previous calculation [JO] in a consistent 
way: we presents results of a calculation of the nuclear Coulomb photoproduction of 
neutral pion pairs, where we include in the elementary amplitude 77* — > 7r 7r° both 
off-shell corrections for the Coulomb virtual photon 7* and contributions from ChPT 
up to 0(p 6 ), and where, in addition, nuclear structure and absorption are taken into 
account in the form factor for the nucleus. 

2. Coulomb photoproduction 

The photoproduction of a 7r°7r° pair in the Coulomb field of a nucleus is schematically 
described by the diagram in figure 1. The virtual photon 7*((/2) for the interaction 
of the incident real photon with the stationary Coulomb field of the nucleus has 
zero energy and transfers only momentum: q 2 = (0, qg). Then the amplitude of the 
reaction j(qi)A — > ir (pi)ir (p2)A has the form 

Tc = 2M A ^F A (q 2t , q^T^^ , (1) 
1 02 1 

where M A and Z are, respectively, the mass and charge of the nucleus. In flU), 
the nuclear form factor F A (q2t, Q2\), which includes nuclear absorption, depends on 
transverse and longitudinal components of the momentum transfer q 2 measured 
relative to the momentum qi of the incident photon (Fa is normalized to Fa{0, 0) = 
1). X^q 7 17 -* is the tensor component of the amplitude of the process 7(q , i)7*(<?2) — > 

From Lorentz and gauge invariances the general parameterization for the 
amplitude T^q 7 _Mrir ' at 0(p 6 ) has the form 

jtrf-Hfi**) = A(s? t} u . g 2 } ^ _ q2)xqiv 

+ B(s, t, u; ql) 25A tl A u - v 2 g^ v - 2u(A^q lu - q 2 ^A u ) 
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+ D(s,t,u; ql) vq 2il q 2v + sA^q 2u - ql{vg^v + 2A fl q lv ) , (2) 

where = (pi-p 2 ) M , s = (gi+g 2 ) 2 = (Pi+P 2 ) 2 , § = s-q%, t = (pi-qi) 2 = (g 2 -p 2 ) 2 , 
u = (P2 - <?i) 2 = (92 - Pi) 2 and v = t — u. 

The differential cross section for the photoproduction of a neutral pion pair in 
the Coulomb field of a nucleus is defined as 

dac = 4eM A (2nf8 El E 2 M A 2 |Tc| d * d P* d * 

where and pj (i — 1, 2) and £ are the energies and momenta of the pions and the 
energy of incident real photon, respectively. For large energies of the incident real 
photon and for small momentum transfers |q 2 | to the recoil nucleus and neglecting the 
offshellness of the Coulomb photon and nuclear corrections, the method of equivalent 
photons allows us to relate the differential cross section for photoproduction of 
pion pairs on nuclei to the total cross section for the process 77 — > tttt: 

^^l = ^ zHo j^Ly_^, , w (3) 

as 7r \2m 7r / s 

where s = vn 2 ^. In this limit equation (^[) enables us to extract model- independent 
information on the process 77 — > tttt from the experimental data on the nuclear 
Coulomb photoproduction of pion pairs. For a more general kinematics, the nuclear 
form factor and the offshellness of the Coulomb photon have to be taken into account, 
however. 

The nuclear form factor i*A(g2t) 921) in © can be estimated in the same 
approximation as in paper since the amplitude of Coulomb photoproduction 



on a single proton at small g 2t is - upon evaluating (Q) - proportional to g 2t : 



In the approach of paper jnj the nucleus is treated as a completely absorbing sphere, 
with the form factor 

F A (qu qi) = qiRMqtR)K x { qi R) + ^^-J^q^K^R) + AF A (q t , q t ) , (4) 

qtR 

where R is the radius of the nucleus, and J n and K n (n = 0, 1) are Bessel functions. 
In our estimates R is chosen to be R = 1.12A 1//3 fm, where A is atomic weight 
of a nucleus. The first two terms in (^) arise from the integration over the three- 
dimensional space outside a cylinder of radius R of the nucleus. They reflect the 
assumption, that the nucleus is completely 'black' for the outgoing pions for impact 
parameters b < R, resulting in a cut 9(b — R), in the profile function. This drastic 
assumption is mediated by the correction term AF A (q t , q{), which arises from the 
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integration over the three-dimensional cylinder behind the nucleus and corresponds 
to the interaction of photons with the nuclear Coulomb field after passing through 
the nucleus. 

In general, the integrations above can be done only numerically. However, when 
both (qiR) and (q t R) are small compared to unity, the correction AF^(q t , q{) can 
conveniently be expanded to obtain 

AF A (q uqi ) = Hi - qi RJ Q {q t R)K x { qi R) - ^^J^q^K^R) 
4 I q t R 

+ \{ qi Rf + [q t Rf] ( - I + [ - qi R + ^(qtR) 2 
3. Chiral lagrangians 

The calculation of the amplitude of 77* — > 7r°7r° at 0(p 6 ) of the momentum expansion 
of ChPT involves tree-level, one-loop and two-loop diagrams of a chiral effective 
lagrangian. The effective meson lagrangian at 0(p 2 ) can be written in the nonlinear 
parameterization of chiral SU(2) x SU(2) symmetry as 

4? = tr ipjr&rf) + ^ tr (tf/t + U X ] ) , (5) 

where 

u{x) = exp G^ (x) ) ' y,= (X- ^° + ) (6) 

represents the pseudoscalar degrees of freedom, F is the bare 7r decay constant, and 
X = diag(x^,x^) is the mass matrix. The covariant derivatives and contain 
the vector and axial- vector degrees of freedom, and are defined as 

D,U = dpU + {Aft - UA*) , = + (AjC/t - E/UJ) , 

where A^ L = ± A^ are the right/left combinations of vector and axial-vector 
fields. The interaction with the electromagnetic field can be included by replacing 
Vn — > Vfj, + ieA^Q, where Q is the matrix of quark electric charges. 

The effective meson Lagrangian of 0(p A ) is presented in the general form with 
structure coefficients Li and Hi introduced by Gasser and Leutwyler in paper @], 

42 = fa - \L 2 ) {txL^f + L 2 tr(±[L„ L u } 2 + 3(L M Z/) 2 ) 
+ L 3 tr((L^) 2 ) - LML,L») tr( X rf + U X ] ) 
- L 5 tr[L^( X tf + U X ] )] + L e (tr( X tf + U X ] )f 



+ L 7 (tr( X f/t - U X ] )f + L 8 tx( X U^ + U X ] U X ] ) 

- L 9 tT(F*R»R" + F^U) - L 10 tr(UF*U^F L ^) 

~ H.tr^Fp 2 + (Fp 2 ) + HMXX}) ■ 

In equation <g), = DJJ-W, R» = W D^U, and F$ L 
right/left combinations of the field strength tensors 

Zv = d,V u - d u V, + [V„ V v ) + [Ap, A v ] , = d„A v - d u A^ + [V„ A v \ + [A„ V v 



(7) 



FV ± are the 



F v 
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It is convenient to present the part of the effective lagrangian at 0(p e ) for 
— > 7i°tt° with structure coefficients di, 



C *-FZ 



di^a^tr (d a U Q d p UlQ 2 ) + d 2 T^ v tx (d a U Q d a UlQ 2 

+ d 3 F^mr (x(U + Ul)Q 2 ) + d^^ti{Q 2 )tv (x(U + tf ) 
+ d 5 ^ a ^ti (Q 2 ) tr (dPUadpUS) + d^^tr (Q 2 ) tr (d a U d a U j ) 
+ diT m T»hx (&*U Q UlQ) tr (dpU UiQ) 

+d 8 F^tr (d a U Q UlQ) tr (d a U tf Q)] , (8) 

where = d^A u — d u A^ is the standard electromagnetic field strength tensor, 
and Uq = exp(i(po/Fo), <po = diag(7r°, —it ). The lagrangian of equation (g) can be 
obtained from the most general representation of the full lagrangian of paper |L7| . 

The structure coefficients of the 0(p A ) and 0(p 6 ) chiral lagrangians of equations 
(0) and (El) can be fixed either from low-energy meson phenomenology, as in papers 
(3J, ffl, or from the bosonization of NJL-type effective quark models (see papers 
T5|, [TBI, ITTJl ITTJ and references therein). 



4. Amplitudes of 77* 



In general, the prediction for the Born amplitude of 77* — > 7r°7r° at 0(p 6 ) involves 
eight structure coefficients di of the lagrangian (|8|). In the NJL model only the 
structure constants di, d 2 , d 3 get nonzero values and contribute to the amplitude: 



^B( P 6) 
B B(p 6 ) 



9F 4 



10e' 
9^ 



di . 



The one-loop amplitude of 77* 
is given as 

^u(p 4 ) = 2e2 1 



0^0 



TV TV 



at 0(p 4 ), which has no UV divergences, 



3F 2 16vr 2 



(Qs-Smi + xt + xi) 



2q 

s 



R 11(p4) = D 11(p4) 
with 

and 



H. 



tt{s, q 2/ 



2m 



l + ^(j;\s)-J-\ q: 



J2(a) 



d x x n log 



ml — ail 



x)x 
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mi 



In our approach UV divergences, resulting from meson loops at 0(p 6 ), are 
separated using the superpropagator regularization method fl8| , which is particularly 
well-suited for the treatment of loops in nonlinear chiral theories. The result 
is equivalent to the dimensional regularization technique used in paper ||, the 
difference being that the scale parameter \x is no longer arbitrary but fixed by the 
inherent scale p, of the chiral theory, namely, jx = Air F . For comparison of these two 
methods, the constants from the UV divergences in the dimensional regularization 
have to be replaced by a finite term using the substitution 



(C-l/e) 



G 



SP 



2C + 1 + - 



-(i og r- 2 (2^ + 2; 



+(3tt = -i+AC+f3iT. 



2=0 



where C = 0.577 is Euler's constant, e = (4 — D)/2, and (3 is an arbitrary 
constant resulting from the representation of the superpropagator as an integral 
of the Sommerfeld- Watson type. The splitting of the decay constants F n and F K is 
used at 0(p 4 ) to fix C SP « 3.0 for F = 90 MeV. 

The one-loop diagrams at 0(p 6 ) involve the structure coefficients Li of the 
lagrangian (|7|) and give the following contributions 

3 2 1 



3F 4 16vr 2 



2L 9 q 2 2 (6s-8mi + xi + x2) 



+ 12L 2 (3s 2 - 8m 2 (s - m 2 )) + 24(2^ - L 2 + L 3 )(s - 2m 
+ 16L 4 (3(s - 2m 2 )( X 2 + X 2 ) - (3s - 4m 2 )tr X 



-16L 5 ml(x 2 u + X 2 d ) 



- 2 ^{jl{s)-Jl(ql))+H^ q l) 



+ 4(L 9 + L 10 ) (6s - 8m 2 n + X 2 U + X 2 d )G„(s) 



+ 12L 2 



s - 2ml - ^(3g 2 4 - 48g 2 m 2 + 32m*) 
\- 2 {ml< - 21g 2 X + 2q\ + 6v 2 (q 2 2 + 12m 2 ; 



25 

- ^<£(<£ + 19m 2 ) - %v 2 q\{ 



£ + 26m 2 ; 



4mt-g 2 4 + |(g 2 2 



10m" 



4e 2 



+ j 2 (ml< ~ Uqtrnl + q% + 3v 2 {q\ + Ami)) 
+ + 14m 2 ) + ^v 2 q 4 2 (q 2 2 + 26m 2 ) 

5 S 

- 2H„(s, q 2 2 ) (sml - Am 4 n + q\m\ - v 2 - ^ 2 (m' + 2q 2 ) 

-^u 2 q 2 2 (ml + q 2 2 ) 
~ ^(s - Ami) + jz(Q2 - W 2 ml - hv 2 ) - ^v 2 q 2 2 
~ ly^t + \c,{s - Ami) 



F 4 16tt 2 ' 
8e 2 



l-|(5 m |-5 2 2 )-|(10m;- 9 l) 

4 



2(Aml-q 2 ) + ^-(10ml-q 2 2 ) 
4 



3F 4 16vr 2 s 3 



L2^\ q 2 2 Jl(s) s 2 + 2§(13m 2 + q 2 ) + (26m 2 + q 2 ) 



+ Jl(q 2 2 ) [s 2 (Aml - q\) - 2~sq 2 2 (8ml + q\) - ql(2Qml + q\) 
- 3q 2 2 H n {s, q 2 2 )~s 2 (s + m%) + ^g 2 2 §(9s + g 2 )} , 

where G w (s) = <7 W + 2J 4 (s) and C w = C SP + log 

Only two-loop diagrams which are factorizable can be calculated analytically. 
Their contribution can be presented in the form 



2e 2 



3F 4 (16tt 



2^2 



Jim%JM) + ^H^q 2 2 )\- 



32 4 _ 3 



+ 2s\Aml - q\) - -s 3 q 2 2 (32ml - 9q 2 2 ) + 2s 2 ql(8ml - 3g 2 ) 
- 2sqi(32mi - qt) - \rnlq* + i§ 3 (4m 2 - q\){ X 2 u + X 2 d ) 



-^JM) 



Qs\Aml -q 2 2 )~ s 4 (20m 4 + Ihmlq 2 - 29g 4 ) 



16m 2 J 4 - ^s 3 q 2 2 (16ml + 166m 2 g 2 + 9q 4 2 ) 



+ 2s 2 g 4 (60m 4 + 91m 2 g 2 2 + 15g 4 ) 
- 2sql(17Qml + AQmlq 2 2 + q A 2 ) + ^mlq 8 2 (27Qml + 35g 2 ) 

+ \ ( 3(5 4 + s')(Aml - q 2 2 ) - 3 S 3 g 2 (36m 2 - 7q 2 2 ) 



24s 2 g 4 (6m 2 - g 2 ) + 6sg 2 b (20m 2 - 3g 2 



+ q Z(36ml-5q 2 2 ))(xl + X 2 d ) 



+ Jl{s){^ 2 Jl{s)-\H^ql))^ 

- 4s 4 (8m 2 + 9g 2 ) + ^s 3 g 2 (32m 2 + 9q 2 2 ) - 12s 2 ql(8ml + q\) 



— m 4 s 3 + 12s 5 
3 



+ 32sm 4 r g 2 6 + ^s~ 3 (3s - 4m 2 )( X 2 + xl) + ^~ 3 (x* + xlf 

6s 5 ml - 2s 4 (20m 4 + 88m 2 g 2 + 19g 4 ) + 32m 4 § 4 
502 



+ 2s 3 g 2 (l36mt-^m 2 g 2 -5g 4 ) 



3 

6 S 2 g 2 4 (96mt + ^m 2 g 2 2 + g 2 4 
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+ 2sg 2 6 (248m 4 r + 43m 2 g 2 + g 4 ) - -m 2 g 8 (57m 2 + q 2 2 ) 

+ ^ ( - 72s 5 + 36s 4 (3m 2 + 8g 2 ) + 245 4 (3s - 4m 2 ) 

- s 3 g^(420m 2 + 433g^) + 3s 2 g 4 (204m 2 + 97g^) 

- 3sg 2 6 (132m 2 + 25g 2 ) + g 8 (96m 2 + g 2 )) ( X 2 + X 2 ) 



+ ^JMUKs) \ - 12s 5 + 2s 4 (12m 2 + 19g 2 ) - ^m 4 s~ 



ys 3 g 2 (16m 2 + 3g 2 ) + 2s 2 g 4 (40m 2 + 9g 



- 2sg 2 6 (16m 2 + g 2 ) + ^m 2 g 8 + 1 (^20m 2 5 3 - 48s 4 + 145s 3 g 2 



7. 



2\2 



- 147s^ 4 + Slsg* - g 8 ) (xl + X 2 d ) - -s~ 3 ( X 2 + X 2 ) 



- ^H n (s,q 2 2 ) 



2 

3' 



--^s 4 (36m 2 +g 2 ) + ^s 3 ^(H2-/ ;; :; -9g 2 > 



+ ^m^5 3 - 2s 2 g 4 (38m 2 + g 2 ) + jjsg 2 6 (40m 2 + gf) 
~ \<<& + ^ 3 (12m 2 - g 2 )( X 2 + X 2 ) 
^ (s 4 g 2 (109m 2 + 3g 2 ) - 48m 4 § 4 + 2m 2 s 3 g 4 - 4m 2 sg 8 - 3m 2 g 2 10 



6s 5 + s 4 (8m 2 - 17g|) 



is 3 g 2 (76m 2 - 45g 2 ) - ^-m 4 J 3 + s 2 g 4 (28m 2 - g 2 ) 
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- sq 6 2 (12ml + q\) + ^m 2 g 2 8 - i(72m 2 S 3 - 48s 4 + 143s 3 g 2 

3 o 

- 141 S 2 g 4 + 45sq 6 2 + q s 2 )( X l + Xd) + ^(xl + xlf 



+4a 

s 5 



23 1 
s 5 mt - -s A mlq 2 2 + 4m 4 5 4 + -s 3 g 4 (3m 2 + 16g 2 ) 



s 2 g 2 6 (20m 2 + q\) + ^g 2 8 m 2 + \m\qf 



3 ^ " 6 



+ 3^ 5 4 (x' + Xd) 



R 21 (p 6 ) = n2i(p 6 ) = Q 

It should be stressed that additional two-loop diagrams, such as box diagrams 
and acnode graphs, which cannot be evaluated analytically, can be neglected: the 
numerical estimates in paper || indicate the smallness of their contributions in the 
photoproduction process under consideration. 

5. Numerical estimates and conclusion 

Monte Carlo techniques were used to compute the total cross sections for the 
photoproduction of 7r°7r°-pairs in the Coulomb field of the carbon (Z = 6) and 
silicon nuclei [Z = 14). It is important to note that for a momentum transfer 
cutoff g max = |q2|max «C e, the effective mass of nn system varies in the range 
4m 2 < ml n < 2e:g max . 

The dependence of the total cross section of the reaction 'yA — > A on the 
momentum transfer cutoff q max for different energies e of the incident real photon is 
shown in figures 2 and 3. (In the calculation, the additional cutoff < 700 MeV 
was used, corresponding to the range of validity of the chiral theory.) In figure 2 
we demonstrate explicitly the influence of the nuclear form factor: our numerical 
estimates indicate that the nuclear effects are mainly saturated by the contribution 
of the two first terms of equation while the correction AFa does not exceed 
10%. For our numerical estimates we used for the parameters Li and di the 
values from (14) and Table 1 of paper [llj obtained from the NJL model. These 



values include the resonance-exchange contributions effectively taken into account 
by integrating out the vector, axial-vector and scalar degrees of freedom. Moreover, 
for numerical comparison with the phenomenological approach to 77 — > 7r°7r°, we 
used the structure constants Lj and d( corresponding to Tables 1 and 2 of paper ||. 
The results of our calculations with the parameters of paper are also shown in 
figure 2. 

The main background for the Coulomb photoproduction of 7r°7r° pairs is the 
double pion photoproduction reaction on the nucleon 7^ — > -k°tt°N via the baryon 
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resonance exchange mechanism [19, 20]. In the photon energy region near e — 1 
GeV the total cross sections of the reactions ^p — > 7r°7r°p and — > 7r°7r°n are 
dominated by diagrams in figure 4 with A(1232)-isobar and resonance iV*(1520) [ [H|| . 
Kroll-Ruderman diagram dominating in other isospin channels and meson exchange 
diagrams do not contribute to the photoproduction of 7r°7r° pairs. The contributions 
of diagrams j, 1, and n of figure 4 to the reaction 7 A — > tt 7t°A can be canceled by 
choosing an isoscalar target because of opposite signs of the corresponding amplitudes 
on the proton and neutron. The amplitude of diagram p in figure 4 is proportional to 
the difference of the coupling constants (gO? — g^) (see paper jn|] for more details). 
The corresponding contribution to the process jA — > 7t°tt A on the isoscalar target 
proves to be also suppressed due to the fact that 

J 0.157 , for proton , 
I —0.136 , for neutron. 



(5? -5? 



The effect of such self-cancellation of diagrams j, 1, n and p from figure 4 is 
demonstrated in figure 5 where we present the results of our calculation of the total 
cross section of the double pion photoproduction reaction on the proton, neutron 
and on the system (proton + neutron). Numerically our calculations are in a good 
agreement with the results of paper |2l| and in a reasonable agreement with the 



recent experimental data [p2 |. 

In figure 6 we also show the photon energy dependence of the total cross sections 
for the coherent double pion photoproduction on the carbon and silicon nuclei 
(isoscalar targets) calculated with the Born nuclear electromagnetic form factor with 
symmetrized Fermi-density [ p3fl : 

. 47r 2 bRpW ( TTfo . „ 

-T d [a) = : — — cos git — sin git cot nbq 

qsmhnbq V it 

where g = |qi — pi — p2 1 - The normalization constant is defined as 

p (o) _ A 1 



4vr R(R 2 + vr 2 6 2 ) ' 

and the following values of radius of the nucleus R and impact parameter b were 
used: R = 2.214 fm, b = 0.488 fm — for carbon, and R = 3.085 fm, b = 0.563 fm - 
for silicon. Our results for cr to4 (C 12 ) agree with recent estimates [24]. 

To estimate the contribution of the reaction 7iV — > 7i°7i°N to this part of the 
phase space which corresponds to the cutoffs introduced above for the Coulomb 
photoproduction of 7r°7r° pairs, the additional cutoffs must be also taken into account 
when Monte Carlo simulating the background: m n7T < 700 MeV and y/—tN < Qmax, 
where = (pi — P2) 2 , Pi and P2 are the momenta of nucleon in the initial and 
final states, respectively. The photon energy dependence of the cross sections for 
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the carbon and silicon targets calculated with the cutoffs are shown in figure 6. 
Our estimates of the 7A — > 7i°7i°N background have shown that in the case of 
an isoscalar target it is strongly suppressed in the energy and momentum transfer 
regions under discussion. According to the symmetry properties of the background 



coherent reaction — > 7r 7r°v4 on isoscalar nuclei p3, the two 7r u mesons prefer 



to propagate together in the CMS. The Monte-Carlo simulation shows that the 
additional suppression of the backgraund to the level lower than the Coulomb 
photoproduction can be easily achieved using some angle cutoffs. To extract the 
signal of Coulomb photoproduction is possible in the experiment where a silicon 
detector is used as a sensitive target to develop a trigger on the nucleus recoil 
momentum. 

Summarizing, we have found that for energies of the incident photon e < 4 GeV 
and for momentum transfers | CJ2I < 200 MeV, 7r°7r° pairs with 2m- K < m niT < 700MeV 
in the reaction 'jA — > 7T°7r°A are produced with a cross section of typically a ~ 40 pb 
and a ~ 120 pb for carbon and silicon nuclei, respectively. Our results demonstrate 
that this reaction can be experimentally investigated with presently available photon 
beams as a new source of the low-energy data on the process 77 — > 7r°7r°. 
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Figure 1. Photoproduction of pion pairs in the Coulomb field of a nucleus. 
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Figure 2. Dependence of the — > vr vr A cross section on the momentum transfer cutoff 
<7max for the energy e = 4 GeV of the incident photon. Compared are: Lj, di from the 
bosonization of the NJL model (JLT]], full curves) and from phenomenology ([H], dash-dotted 
curves), without and with a nuclear form factor (indicated by (1) and (2)). 
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Figure 5. Total cross sections of the reactions jp — > 7r°7r°p (full curve) and jn — ► 7r°7r°n 
(dashed curve) calculated in the framework of baryon resonance model [^] with diagrams 
of figure 4. Dash-dotted curve corresponds to the total cross section for the sum of 
the amplitudes on the proton and neutron. The experimental points for the reaction 
7P — ► tt°tt°p were measured at DA<I>NE \22\. 
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Figure 6. Total cross sections of the double pion photoproduction reaction (atot) and 
cross sections calculated with the cutoffs y/—tjsr < 150 MeV and m n7T < 700 MeV (a cu t) 
for coherent processes on the carbon (dashed curve) and silicon nuclei (full curve). 



